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In this paper, phase transitions in the Mitsui model without longitudinal field but with a transverse one are 
investigated in the mean field approximation. The one-to-one correspondence has been established between 
this model and the two-sublattice Ising-type model with longitudinal and transverse fields. Phase diagrams and 
diagrams of existence of the ferroelectric phase are constructed. In the case Si = (SI is the transverse field), 
a simple analytical expression for the tricritical temperature and the condition of existence of the tricritical point 
are obtained. For Si jt 0, systems of equations for the tricritical point and for the condition of its existence are 
written. 
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1. Introduction 

The Mitsui model was proposed in 1958 [1J to theoretically explain the ferroelectric properties of 
the Rochcllc salt. In 1971, Zeks, Shukla, and Blinc [2 formulated this model in terms of pseudospins 
and just in this form it is known at the present time. 

For better quantitative description of the Rochelle salt the conventional Mitsui model was 
modified in many ways. For instance, in reference [3] an additional piezoelectric interaction was 
included in the Hamiltonian of the model and in references [I] and [5j a transverse field (or tun- 
neling) was included as well. To take into account a realistic structure of Rochelle salt crystal, 
the four-sublattice Mitsui model was considered j6j. Besides the Rochelle salt, the Mitsui model 
with transverse field was used for theoretical description of some other ferroelectric compounds, 
notably, RbHS0 4 [H|S] and NH 4 HS0 4 crystals. 

In reference |S] the effect of hydrostatic pressure on thermodynamic properties of the Rochcllc 
salt was studied. Hydrostatic pressure makes it possible to change the parameters of the model. 

It turns out that the Mitsui model covers not only the ferroelectrics of the order-disorder 
type but also other physical objects with two-minimum asymmetric potential. For instance, in 
references [TU] and [TT] pseudospin-electron models based on the Mitsui model were considered. 

Despite a rather wide use of the Mitsui model, there is no detailed analysis of its phase behavior. 
In reference |T5] the diagram of existence of the ferroelectric phase was calculated. However, this 
diagram is far from being complete. More detailed though still incomplete diagram was obtained 
in reference |13| . 

In the present paper, a rigorous and original mathematical investigation of phase transitions in 
the Mitsui model is proposed. We managed to construct a complete phase diagram of the Mitsui 
model in the mean field approximation, first without tunneling (or transverse field) and then with 
nonzero tunneling. For all curves (surfaces) of the diagram the analytical expressions, equations or 
systems of equations are given. 
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2. Hamiltonian of the Mitsui model in the mean field approximation 

The Hamiltonian of the Mitsui model without external field reads 

H = \ E J » & As f + s * Bs f) - E K v% A sf - a E ( S * A - S * B ) - n E ( S ? A + S * B ) ■ 

ij ij i i 

(2.1) 

Here subscripts A and B denote two sublattices, Sf is a-component of the pseudospin on ith site 
of sublattice A; Jij and are the interactions between pseudospins of the same sublattice and of 
different sublattices, respectively, A is the asymmetry of the anharmonic potential; transverse field 
SI describes the tunneling between two wells of the two-minimum potential. The model considered 
is a lattice one. However, since we use the mean field approximation, we do not need to specify the 
type of the lattice. 

The transformations Sf — > —Sf B , Kij — > —Kij transform Hamiltonian (2.1l into the Hamil- 
tonian of the two-sublattice model with longitudinal field A (the same for both sublattices) and 
transverse field Q. Therefore, the two Hamiltonians are equivalent. 

In the mean-field approximation, Hamiltonian (2.1l reads 

H = y U VaVb + ^ + ,D}+E W + H ?) > ( 2 - 2 ) 

where the following notations are introduced: 

Ht = -(A + K Vb + J Va ) Sf A - nSf A , (2.3) 
Hf = (-A + K Va + J Vb ) Sf SlSf, (2.4) 

K = E^«- = E K » > J = E J ^ • = E J ^ > (2- 5 ) 

i j i j 

r\A = (Sf A ) and r\B = (Sf) are the average values of pseudospin on sublattices A and B, respec- 
tively, N is the total number of pseudospins. 

3. n = case 

3.1. Free energy, thermodynamic equilibrium conditions and order parameters 

Let us first consider the SI = case. In this case the eigenvalues of Hamiltonians Hf and Hf 
are as follows: 

K =-r(A + K7? B + Jr, A ), -X A ; 



\ B =--(-A + Kr) A + Jr) B ), -X, 



(3.1) 



The partition function for one unit cell reads 

Zi = (e-^ + eP x *) (e-^ B + AJe^K'^^W)] . (3.2) 
The free energy per one unit cell is as follows: 

F = -9 In Zi , 

F = Kr) A r) B + ~ ( V 2 A + r) 2 B ) - 8ln (e - ^ + / A *) -01n(e^ A « +eP X *), (3.3) 
where 9 = is thermodynamic temperature. From thermodynamic equilibrium conditions 

OF 
9Va 



dF 
dr) B 



= 0, 

T - A (3.4) 

= 

T.A 
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we obtain the following equations for r\ A and r\ E 



2V A 
2 Ve 



tanh (- 
tanh (• 



-/3AJ, 



(3.5) 



Taking into account equations (3.5), we can rewrite the expression for free energy in the following 
form: 



J , 

F = Kr) A r) B + - (t} 



-In 
2 



(3.6) 



The system of equations (3.5 1 has the solutions of two kinds: for the first ones rj A = —rj E 



(then the system is reduced to one equation only), they exist for arbitrary A; for the second ones 
r) A =/= — r\ B , they exist in a bounded region of values of A and correspond to the ferroelectric phase. 

Let us divide the (K, J) plane into eight segments as shown in figure [T] The plots of the free 
energy as a function of A at zero temperature in the centers of the unit-circle arcs for every segment 
are shown in figure [2] (It is easy to obtain these plots from equations ( |3.5| and (3.6) setting 9 — >• 
or f} — > +oo.) As one can see from figure 2, in segments 4 and 5 there are no phase transitions, 
in segments 6, 7, and 8 there are phase transitions at A = only. In segment 3 there are only 
second-order phase transitions that are easy to investigate. The most complicated and interesting 
picture of phase transitions is observed for segments 1 and 2, therefore we consider only the region 
K > 0, J > 0. 




Figure 1 . Division of the (K, J) plane into eight segments, which correspond to different relations 
between parameters K and J. 



Let us introduce new variables: ferroelectric order parameter £ = r\ A + r\ B and antiferroelectric 
order parameter a = r\ A — r\ B . Let us also divide all energetic values by K + J and introduce the 
following notations: 



a 



K-J 
K + J 



7 



K + J 



t : 



K + J 



f 



F 



K + J 



Now equations (3.5 1 can be rewritten in the following form: 



(i + 2 -^ 2 
(i-0 2 -a 2 



3 -2 7 /i 



(i 



(3.7) 



(3. 



If £ = (first type solutions, both sublattices are equivalent up to a sign of 77) then the first 
equation becomes an identity and the system (3.8 I is reduced to a single equation. If £ 7^ (second 



23603-3 



Yu.l. Dublenych 




-0.7 

-0.8 




0.7 — | jp 




0.7 -| F 




0.7 -| F 




0.7 -| F 



-0.8 0.0 





A 

i 

0.8 



-0.7 - 
-0.8 




Figure 2. Free energy as a function of A at zero temperature. Numbers over figures correspond 
to the unit-circle points indicated in figure [l] 



type solutions) then equations (3.8 1 can be rewritten in a simpler form: 



1 + c«/* 
1 - e«/* 



7 = _ (7+ _ +tln __ 



(3.9) 



In view of the symmetry we consider only a 0, £ and 7^0. 
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The expression for the free energy per unit cell in terms of new variables reads 

/ = \ (e - aa 2 ) + * In{ [l - (£ + af] [l - (£ - a) 2 ]) - 2tln2. (3.10) 

3.2. Second-order phase transitions 

The second-order ferroelectric phase transition corresponds to the branchpoint of the curve 
17(7) (a and t being fixed) or of the curve a(t) (a and 7 being fixed) where a solution of the first 
type turns into a solution of the second type. We denote the value of a in this point by o. It is as 
follows: 

a = lim a = Vl - 4f . (3.11) 
One can see from this expression that the second-order phase transitions exist up to the temperature 

^max — ^ • (3.12) 



Substituting £ = and a from equation (3.111 in equation (3.9 I, we obtain the equation for the 



curve 7 (i) of the second-order phase transitions: 

7 =£* + ilni±*. (3.13) 

Z l — O" 

Now let us find the minimal temperature for the existence of the second-order phase transitions 
at fixed a. If a = 1, then only second-order phase transitions exist. If — 1 < a < 1, then there are 
phase transitions of both second and first orders. The latter exist from t = to a certain value of 
the temperature. As one can see from figure [3j the tricritical point, i.e. the point where the order 
of phase transition changes, can be determine from the following condition: 

lim-^ = 0. (3.14) 
do- 




Figure 3. Antiferroelectric order parameter as 
a function of 7 for several values of temper- 
ature (numbers over the curves), a = 0.0875. 
Therrnodynamically stable states are depicted 
by heavy lines. At t — 0.03 and t = 0.1 there are 
first-order phase transitions and at t — 0.1507 
and t — 0.2 there are second-order ones. 



Figure 4. Antiferroelectric order parameter as 
a function of 7 for several values of parame- 
ter a (numbers near the curves). The values 
of the temperature are calculated using expres- 
sion (3.151. At a = —0.1 and a — 0.25 there are 
second-order phase transitions and at a — 0.35 
and a = 0.45 there are first-order ones. 
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We obtain for the tricritical temperature: 



ttc — 



6(a- 1) ' 



(3.15) 



The tricritical point exists if, at the temperature determined by equation (3.15), the following 
condition is satisfied: 

d 2 7 

lim —L < ( 7 > 0). (3.16) 

£-s-o da 1 

This is clear from figure [4] where curves a( r y)/2 for several values of parameter a and corresponding 
values of the temperature [see equation (3.15 1] are depicted. The existence of the region where the 
dependence 17(7) for £ ^ is two-valued indicates that there is a first-order phase transition. This 
two-valuedness disappears with decreasing a and the order of the phase transition changes. From 



equations (3.15) and (3.161 one obtains: 



a ^ 



(3.17) 



equation (3.131 yield 



Let us find the maximum of curve 7(4) at fixed a. The extremum condition d'y/dt = and 

(3.18) 



In 



1 + or 



1 — <J 



a + 1 
a 



+ (a - 1)5- 



Excluding a, we obtain the equation for 7: 

(27 + y/ {2<y)> -a> + i ) tanh WgEg±j = fl + 

3.3. First-order phase transitions 



1. 



(3.19) 



Up to here we analyzed the second-order phase transitions and found the expression for the 
surface 7 = j(t, a) of these transitions as well as the tricritical point and the condition for its 
existence. Now let us consider the first-order phase transitions. 

A first-order phase transition corresponds to a point of self-intersection of the curve for the 
free energy /(£, a) as a function of 7 at fixed a and t (or as a function of t at fixed a and 7). But 
only this point of self-intersection gives the first-order phase transition in which the multivalued 
function for the free energy takes the minimal value from all possible values at fixed 7. 

If —1 ^ a ^ 1/4, then the first-order phase transitions exists for ^ t < t tc [see equation 
(3.15)]. The curve for them in the (7, i)-plain (i.e, the intersection points of the branches £ = 



and £ 7^ of the free energy) can be found from the following system of equations: 



1 + c^ 

1 - ef/' ' 

7 = — it H hi In 

' 2 2 1-cr + C 

a l + o"! 

7 = 7,0-1 +tln- , 

I 1 — (J\ 

f(Z,<r)=f(0,<r 1 ). 
If t = 0, then the first-order phase transition occurs at 

7 = ±~(a + l). 



(3.20) 



(3.21) 



In the t ^ case the system of equations (3.201 can be solved only numerically 



23603-6 



Phase transitions in the Mitsui model 



The phase coexistence curves for several values of parameter a are shown in figure [5] The 
region of ferroelectric phase is bounded by such a curve and by the coordinate axes. The curve of 
tricritical points [more exactly its projection on the plane (7, t)] is also shown in figure [5] (heavy 
line). From this curve another two curves branch off: the curve of minima of function j(t) for all 
possible values of a and the curve of maxima for the second-order phase transitions. The curve of 
tricritical points furcates into the curve of branchpoints and the curve of critical points (dashed 
line). If a > 0.1793293 (four curves on the right), then there is a region with two second-order phase 
transitions (see also figure [61. At this value of a the maximum of the curve j(t) of the second-order 
phase transitions coincides with the tricritical point. 




Figure 5. Coexistence curves in the (7, f)-plane for several values of parameter a (numbers near 
the curves). The curve of tricritical points (heavy line), the curve of minima for the first-order 
phase transitions, the curve of maxima for the second-order phase transitions, the curve of 
branchpoints and the curve of critical points (dashed line) are indicated. 



If a ^ 1 /4, then the upper part of the curve of the first-order phase transitions corresponds to 
the phase transitions within the ferroelectric phase. The temperature of these transitions (at fixed 
a and 7), i.e. the temperature for the self-intersection points of the free energy curve, can be found 
from the following system of equations (the solutions of the type £ = £1 should be rejected): 




l-<7 + £ 



1+efr/t (3.22) 

a £1 . l + <7i -6 
7 = 2 <Tl + T + * ln l-a 1+ 6' 

/&ff)=/Ki,*i). 



In figures [6|b) and (c) some fragments of phase diagrams in coordinates (7, t) for a ^ 1/4 
are shown. The upper part of the curve of the first-order phase transitions corresponds to the 
transitions within the ferroelectric phase. 

If a ^ 1/4, then the critical temperature can be determined from the following system of 
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Figure 6. Fragments of the phase coexistence curves for a) a = 0.25 (the tricritical point is 
shown), b) a = 0.46 and c) a = 0.7 (see figure [H|. 



equations: 



6*=°> 
da 2 ' 



/ 1 + e^/* 

'l + £ 2 + 2£- 



1 - e «/* ' 

7 = — <T H hi In . 

' 2 2 1 -cr + C 

This system of equations can be rewritten in the form: 

1 + £ 2 - 2t + - (t - yja 2 e + 2ai£ 2 (l - a) + t 2 (l + a) 2 ) 



1 + < 2 - 3< + - ( i - J K 2 (l - a) + t(l + a)] 2 - 4a 2 £ 2 (2ai - 2t - a) 



a 2 = 


1 + £ 2 - 2t + - 
a 


a 2 = 


1 + < 2 - 3< + - 


a 2 = 


o 1 + 

l + £ 2 + 2£ 

1 - 


7 = 


a £ , 1 
2 (7+ 2 +<ln I 



whence it follows 



(3.23) 



(3.24) 



-g - [(3a - l)t -a]y/q- 2(a + l)t [(a - l) 2 t - a 2 ] 
5 2o(o-l)[2(o-l)*-a] ' 1 j 



23603-8 



Phase transitions in the Mitsui model 



where q — (a— l)(3a 2 — 14a— l)t 2 — 2a(2a 2 — 5a + l)t + a 2 (a— 1), and hence the system of equations 
is reduced to a single transcendental equation. 

As one can see from figure [6j the curve 7(4) of the first-order phase transitions has a minimum 
if a is big enough. To find it, let us differentiate the last equation of system (3.20 1 with respect to 



t and then substitute dj/dt — 0. After simple transformations we obtain the following equation: 



£ 2 - aa 2 + A-fa 



-aa x + 4701 , 



(3.26) 



of equations (3.201 also satisfy it 



which together with equations (3.20 I gives the point of minimum. At zero temperature the solutions 



hm ^ = 0. 

t-*o dt 



(3.27) 



The branchpoint of the curve (at fixed a ^ 1/4) can be found from equations (3.201, substituting 

1-af 



t = 



(3.28) 



3.4. Regions of existence of the ferroelectric phase 

To conclude, we obtained explicit expressions or system of equations for all special points of 
the phase curve at fixed a. This makes it possible to construct a diagram of the regions where the 
ferroelectric phase exists. There are seven regions in the (7, a) plane (figures and D . They are 
bounded by the following curves: 

(1) the curve of minima for the first-order phase transitions [equation ( 3.26[ )j: 

(2) the curve of tricritical points [equations (3.11), (3.131, and (3.15 1]; 



(3) the curve of branchpoints [equations (3.20 

(4) the curve of critical points [equation ( 3.24 1] 



and (3.281] 



(5) the curve of maxima for the second-order phase transitions [equation (3.19 1] 



(6) the straight line of the first-order phase transitions at zero temperature [equation (3.21 1] . 

In the point (a s» 0.179329, 7 = 0.283995) the curve of minima for the first-order phase 
transitions and the curve of maxima for the second-order phase transitions branch off from the 



0.36 -i a 




0.0 0.1 



0.2 0.3 0.4 0.5 0.6 



Figure 7. Regions of existence of the ferroelec- 
tric phase (see also figures [9] and [lO]). Dashed 
line (prolongation of the line which bounds re- 
gion V) is shown only for comparison with the 
diagram from |12j . 



Figure 8. Regions of existence of the ferroelec- 
tric phase (fragment). The filled circle is the 
branchpoint. 
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curve of tricritical points and in the point (a = 0.25, 7 = 0.307041) the curve of critical points 
branches out into the curve of branchpoints and the curve of critical points. All curves except for 
the curve of maxima for the second-order phase transitions converge at a = 1. 



a/2 



f/2 






a/2 



f/2 




Figure 9. Parameters of ferroelectric and antiferroelectric ordering for different regions where 
the ferroelectric phase exists. Only thermodynamically stable states are depicted, a — 0.36. 
From top to bottom: I) 7 = 0.337; II) 7 = 0.3375; III) 7 = 0.338; and IV) 7 = 0.339. For region 
III the phase transition within the ferroelectric phase is indicated by filled circles. 



In figures [9] and 10 the behavior of parameters <r and £ is shown for every region. In region 
I (figures [7] and [sj) only one high-temperature second-order phase transition exists. In region II, 
in addition to the mentioned one, there are two first-order phase transitions. In narrow region III 
there are four phase transition: two of the first and two of the second order; one first-order phase 
transition is within the ferroelectric phase. In region IV there are three phase transitions, one of 
them being the first-order transition. In region V there are two second-order phase transitions just 
like for the Rochelle salt. In region VI there is only one phase transition (of the first order) and in 
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Figure 10. (Continuation of figure [9]) V) a = 0.36, 7 = 0.34; VI) a = -0.4, 7 = 0.14; 
VII) a = 0.36, 7 = 0.35. 



region VII there are no phase transitions at all. 



4. Q ^ case 

4.1. Free energy and conditions for thermodynamic equilibrium 

Now let us consider the case of nonzero tunneling. The eigenvalues of one-site Hamiltonians 



(2.3 1 and (2.4 1 are as follows: 

1 



-K , A a - A + K Vb + J Va ; 



A, -ZJAI-VJ 

A B =-^/A; -X,. A„ = -A +/v n t + ./,/. 



(4.1) 



and the free energy per one unit cell reads 

f = Kt, aVb + ^(vl + nl) - \ (v 7 ^ + ™ + 

- 01n(l + e-' 3 A^ W )-01n(l+c^V S J W ). 



S2- 



(4.2) 
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The conditions of thermodynamic equilibrium (3.4 1 yield the following equations: 

A, 



2q 



tanh (— j8A A ), 



2ri R = . Ab tanh(-/3A„). 

Like in the SI = case, let us pass to dimensionless values, introducing one more notation: 





(4.3) 



Let us also introduce new variables: 

A. +A, 

x 



y = 



2{K + J) 
A A -A B 
2(K + J) 



K + J 



2 2 



= 7 — a- 



(4.4) 



(4.5) 



= 7 — a- 



Then the system of equations (4.3) becomes: 



2x = A(x + y)~ A(-x + y), 
2(7-2/) 



A{x + y)+A{-x + y), 



where 



A(z) 



2Vz 2 + uj 2 



tanh 



\/z 2 + UJ 2 
2lT 



(4.6) 



(4.7) 



and the expression for free energy takes the form: 

(7 - y? i 



/ = x 2 - ^— VL - - (v^ + y) 2 + cu 2 + y/(x - y) 2 + lu 2 ) 
- t In (l + c -V(x+v) 2 +" 2 /^ - t In (l + e - V(*-«) 2 +" a /*) . 



(4.8) 



4.2. Second-order phase transitions 



Like in the f2 = case, the system of equations (4.6 1 has solutions of two types: 1) x = 0; then 
the first equation becomes an identity; and 2) x ^ 0, which corresponds to the ferroelectric phase. 
Letting x tend to zero, we obtain from (4.6 1 the system of equations for hypersurface 7 = 7(0;, a, t) 
of the second-order phase transitions: 



(4.9) 



B{y) -1 = 0, 
7 = y + aA(y), 

where B(z) = d^4(z)/dz, y — lim y. Having determined B(z), we can write the system in the form: 



2tuj 2 s 



y 2 Vy 2 + uj2 



At (y 2 + uj 2 ) i 



y 2 



0. 



7 = y + 



asy 



(4.10) 



2^jpT^ 
where the following notation is introduced: 

s = tanh 



Vy 2 +< 



2t 
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Setting 7 equal to zero in equation (4.101, we obtain the expression for maximal temperature at 
which the phase transitions exist: 

imax = 1 + ■ ( 4 - U ) 
In 



1- 2w 

The logarithm in the latter expression makes sense if uj ^ 1/2. Hence, 

1 



(4.12) 



Differentiating equations (4.10) with respect to t, setting dj/dt equal to zero, and eliminating 
dy/dt, we obtain after simple transformations: 

2(1 + a) (f + a, 2 ) (y 2 Sv /^W - t (f + c 2 )) 

= a (ztuj 2 s^y 2 + uS 1 - 2y 2 s^y 2 + uj 2 + 6tuj 2 (y 2 + uj 2 ) - uj 2 y 2 s 2 ^) . (4.13) 

This equation together with equations (4.10) determines the point of maximum for the curve j(t) 
of the second-order phase transitions (at fixed ui and a). 

The equation for tricritical point is similar to that in the = case: 

lim p. = 0. 

k-»o dy 

From this equation we obtain: 

(a+l)D(y)-3a [C (y)] 2 = 0, 
where the following notations are introduced: 

d 3 A(z) 



D(Z)= , , 

dz 6 



The tricritical point exists if both equation (4.141 and the following condition are satisfied: 

d 2 7 

lim —L sC ( 7 > 0). 
x^o dy 2 

This yields the equation 



8 [D (y)} 3 - 9C (y) D (y) E (y) + ? [C (y)f F (y) = 0, 

5 



where 



E(z) 



d 4 A(z) 



F(z) 



d 5 A(z) 



(4.14) 
(4.15) 
(4.16) 

(4.17) 

(4.18) 
(4.19) 



dz 4 ' w dz 5 ' 

If the second-order phase transitions exist until t — 0, then it follows from equations ( |4.10[ ), 
that for this zero-temperature transition 7 is the following: 



7=(2c)- 2 / 3 (a + (2 W ) 2 / 3 ) y, 



where 



^(2 W ) 2 / 3 (l-(2 W ) 2 / 3 ) 



1/2 



Rewriting equation (4.15) for zero temperature we obtain: 

_ 5(2w) 2 / 3 - 4 



(4.20) 
(4.21) 

(4.22) 



4(2w) 2 / 3 - 5 ' 

If a is bigger than this value (at fixed ui), then the second-order phase transitions begin at zero 
temperature. It follows from equation (4.18) that equation (4.221 is satisfied under condition 

uj ^ 2~ 5/2 « 0.1768; 

then a < 1/2. 



(4.23) 
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4.3. First-order phase transitions 

The points of the first-order ferroelectric phase transitions can be calculated from the following 
system of equations: 

2x = A(x + y) - A(-x + y), 
= A(x + y)+A(-x + y), 
A(yi), 



a 

7-yi 



f(x,y) = f(0, Vl ). 

Differentiating the last equation with respect to t and eliminating the derivatives, we obtain the 
following equation: 



{i-vf 



1-V\ (x + y) 2 +uj 2 



x + y 



7 - y\ (x - y) 2 + uj 2 



(7-yi) 2 2( 1 -y 1 )yl+u? 



x-y 



yi 



(4.25) 



which, together with the system of equations (4.241, yields the extremum of the curve ^(t) of the 
first-order phase transitions (at fixed uj) . If t — 0, then the solutions of the system of equations ( 4.24 1 
satisfy equation (4.25 1 as well. Hence, for the curve j(t) of the first-order phase transitions we have 



hm p.=Q. 
t-X) dt 



(4.26) 



Like in the f2 = case, we can obtain the system of equations for the tricritical point setting 
d^/dy and d 2 7/dy 2 equal to zero: 



B_ - Ba 



2B-B+ - B- 



B A 



C+ (1 - B_) 3 + C_ (1 - B + ) 3 = 0, 



2x = A + - A_ 



(4.27) 



where A± = A(±x +y), B± = B(±x + y). a nd C± = C(±x + y) 
tern 

( 2 ( 



Rewriting the system of equations (4.271 for zero temperature, we obtain: 

3 



2 +w 2N3/2_ w 2 



(2 (z 2 + 



9\3/2 9 

w — or 



{zl+uj 2 )' 



2Jz 2 +uj 2 



z++Z-= 0, 



2Jz 



4(z 2 + w 2 ) 5 (z 2 _ + uj 2 )' 2 -uj 4 
co 2 (uj 2 - (4 +uj 2 Y - (z 2 _ + uj 2 ) 



(4.28) 




where z± = ±x + y. The nontrivial solution of this system at fixed uj ^ 2~ 5 / 2 corresponds to the 
point where the curve of the first-order zero temperature phase transitions and the curve of the 
critical points end. 
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Figure 11. Phase coexistence curves in (7, t)-plane for several values of a (number near curves), 
curve of tricritical points (hard line), curve of extrema for the first-order phase transitions (not 
hole), curve of maxima for the second-order phase transitions, curve of branchpoints and curve 
of critical points (dashed line), uj = 0.1. 

The phase coexistence curves in the (7, i)-plane for several values of a as well as the curves of 
tricritical, critical, and branchpoints at ui = 0.1 are depicted in figure [XT] The latter two curves 
do not converge at zero temperature and, hence, there is an interval of values of a where the 
phase diagrams look as shown in figure [12] (b) , i.e., they are composed of the curve of second-order 
phase transitions and of the curve of first order phase transitions within ferroelectric phase. At 
ui = 2~ 5 / 2 the curves of critical points and branch- points converge again at zero temperature, and 
if 2" 5 / 2 < ui < 0.196815 these curves converge at nonzero temperatures and pass into the curve 



a) 




0.015 - 



0.010- ; ; 

0.005 - ; ; 

0.000 -I 1 r 1 -^ \ 1 1 1 -1 

0.400 0.405 0.410 0.415 0.420 

Figure 12. Fragments of the phase coexistence curves at lj = 0.1 (see figure [TT| . a) a — 0.34 (cir- 
cles are the extremum points), b) a = 0.662 (dashed line on the left is the curve of branchpoints 
and the line on the right is the curve of critical points). 
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of tricritical points (figure 13). The latter is many- valued in low-temperature region. This leads to 
the existence of a new type of phase diagram with two tricritical points (figure 14). If uj is bigger 
than ?s 0.196815, only the curve of tricritical points remains. 




0.6 



Figure 13. Phase coexistence curves in (7,t)-plane for several values of a (number near curves). 
uj — 0.18. The heavy line is the line of tricritical points. The curve of tricritical points and the 
curve of branchpoints are indistinguishable on this scale. 



0.0250 



0.0125 



0.0000 



0.35137 



0.35147 



0.35157 



Figure 14. Fragment of the phase coexistence curve at a 
Dashed line is the curve of tricritical points. 



0.4941, uj = 0.18 (see figure 13 1 



4.4. Regions of existence of the ferroelectric phase 



In figures [I5}j20] the diag ram of the regions where the ferroelectric phase exists is shown for 
uj = 0.1. The diagram is composed of the same curves as in the 51 = case as well as of the 



straight line of second-order zero-temperature phase transitions (4.20) and the curve of maxima 



for the first-order phase transitions (4.251 which passes very closely to the curve of first-order 
zero-temperature phase transitions. 

Like in the uj = case, two curves branch off (in the same point) from the curve of tricritical 
points: the curve of minima for the first-order phase transitions and the curve of maxima for the 
second-order phase transitions. If uj < 0.196815, then the curve of tricritical points bifurcates in 
the curve of critical point and the curve of branchpoints which, if uj > 2~ 5 / 2 , conflow again into the 
curve of tricritical points. The curve of first-order zero-temperature phase transitions is composed 
of two parts. The first one corresponds to the transitions from (or into) the ferroelectric phase. 
It begins at the origin of coordinates and ends at an extremity of the curve of the branchpoints 
(uj ^ 2~ 5 / 2 ) or of the tricritical points (lj > 2 -5 / 2 ) (this extremity is an end point of the curve for 
the second-order zero-temperature phase transitions). The second part exists only at uj < 2 -5 / 2 , 
corresponds to the phase transitions within the ferroelectric phase and ends at an extremity of 
the curve of the critical points (see figure 16). In the uj ^ 2~ 5 / 2 case, the point where the curve 
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Figure 15. Regions of existence of the ferroelec- 
tric phase, u = 0.1. 



Figure 16. Regions of existence of the ferro- 
electric phase (fragment), u) = 0.1. (See also 
figure 21 1. 



of tricritical points and the curve of the first-order zero-temperature phase transitions meet is 
determined by equation (4.221; in this point the latter curve smoothly turns into the curve of the 
second-order zero-temperature phase transitions. 

At nonzero uj new regions occur. For instance, in region VIII there are one low-temperature 
first-order phase transition into the ferroelectric phase and one second-order phase transition. In 
region IX three phase transitions exist: first-order one and two second-order ones (the first-order 
phase transition occurs within the ferroelectric phase). In region X there are one first-order phase 
transition within the ferroelectric phase and one second-order phase transition. 

The curves of minima and maxima for the first-order phase transitions converge in region VIII 




Figure 17. Regions of existence of the ferroelec- 
tric phase (fragment), ui = 0.1. The filled circle 
indicates the point where the curve of tricriti- 
cal points forks into the curve of critical points 
and the curve of branchpoints. 



Figure 18. Regions of existence of the ferroelec- 
tric phase (fragment), uj = 0.1. 
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Figure 1 9. Regions of existence of the ferroelec- Figure 20. Regions of existence of the ferroelec- 

tric phase, cj — 0.18. trie phase (fragment), ui = 0.18. Dotted line is 

the upper part of the curve of tricritical points. 



(a rj 0.3433, 7 « 0.3255, t » 0.0484) cutting off a long spit from it (figure [18}. This is region XI 
with two low-temperature first-order phase transitions and one second-order phase transition. 

Further, the curve of maxima for the first-order phase transitions cuts off narrow strips from 
regions IX, V, and VII. These strips are regions XII, XIII, and XIV, respectively, where, in addition, 
two close low-temperature first-order phase transitions appear. In figure [17] these regions are not 
seen because they are too narrow. 

As one can see, at a sufficiently small value of u> the diagram is richer than at u = but if 
the value of w exceeds some number, the regions of the diagram disappear one after another and 
the diagram becomes poorer. Region III disappears firstly. In figure [17] it looks like a short curve 
segment at the beginning in the branchpoint. It is region XII that disappears the next, and, at 
to = 2~ 5 / 2 , region X becomes a point. 



5. Conclusions 



Hence, we performed a complete analysis of phase transitions in the Mitsui model (without 
and with transverse field f2) in the mean field approximation. Some results concerning the phase 
diagram of the Mitsui model were obtained earlier [5J [TSJ [TB] but they were incomplete (for the 
O = case) or partial (for the il 7^ case). 

In the f2 = case, we derived an analytical expression for tricritical temperature and the 
condition of its existence. In this case, there are seven regions in the plane (a, 7) that correspond 
to seven different types of behavior of order parameters. At sufficiently small but nonzero Q their 
number doubles. With O increasing these regions change their form and shift in the (a, 7)-plane. 
Starting with certain value of ft the number of regions decrease and at u ^ 1/2 only the region 
without phase transitions remains. 

At nonzero f2, second-order phase transitions are possible at zero temperature, which is not 
possible at fl = 0. The maximal number of phase transitions is four and five at ft — and at 
fi / 0, respectively. 
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a/2 



a/2 



0.0 0.1 0.2 

a/2 



0.0 0.1 



a/2 




0.0 0.1 




Figure 21. Ferroelectric and antiferroelectric order parameters for the regions of existence of 
ferroelectric phase, lu = 0.1. Only thermodynamically stable states are depicted. From top to 
bottom: VIII) a = 0.65, 7 = 0.4035; IX) a = 0.66, 7 = 0.4065; X) a = 0.667, 7 = 0.4085; and 
XIV) a = 0.2, 7 = 0.2912765. The first-order phase transition within the ferroelectric phase is 
shown by filled circles. 
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0a3OBi nepexoflu b Mop,en\ Mii4yi 

rO.I. fly6^eHMH 

I HCTMTyT 4>i3WKH KOHfleHCOBa h mx cucTeM HAH YKpaTHM, Byji. I. CBeHL4ii4bKoro, 1, 79011 JlbBiB, yKpama 

B po6o~ri b Ha6/in>KeHHi cepeflHboro nojia flocnifl>KeHO cpa30Bi nepexoflw b MOflejii MiLiyi 6e3 no3flOB>KHboro 
nojia, npcrre 3 nonepeHHHM ncneM. BcTaHOB/ieHO B3aeMooflH03HaHHy 3ajie>KHicTb Mi>K TaKOKi MOfle^jiio Ta 
flBoniflrpaTKOBOio MOfleji/iKD Tuny bmra 3 no3flOB>KHiM i nonepeMHi/iM ncxnnMM. rio6yflOBaHO cpa30Bi fliarpaMH 
Ta fliarpaMU o6jiacTe£n icHyBaHHa cerHe~rocpa3n. fljia BunaflKy Q = (Q — nonepeHHe no/ie) oflepxoHO 

npOCTM M aHajliTHHHHH BMpa3 p,flft TpM Kpi/ITWHHO'f TeMnepaTypi/1 V\ yMOBy icHyBaHHfl Tpi/lKpMTI/lHHOI tohkh. fl/in 

Q ^ 3anncaHO cucTeMH piBHSHb fljia TpuKpnTMHHOi tohkm v\ ywioBH TV icHyBaHHa. 
K/iKJHOBi enema: dpa30Bww nepexifl, cerHerocp33a, MOfle/ib Miijyi, Tpi/iKpi/m/i hhs tohks 
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